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We study the ground-state phase transitions of quasi-one-dimensional quan- 
tum Heisenberg antiferromagnets by the quantum Monte Carlo method with 
the continuous-time loop algorithm and finite-size scaling. For a model which 
consists of S = 1 chains with bond alternation coupled on a square lattice, 
we determine the ground state phase diagram and the universality class of 
the quantum phase transitions. 

1.1 Why Quasi-One-Dimensional Systems? 

Low-dimensional quantum antiferromagnets have attracted much attention 
in recent years. Due to quantum fluctuations, they often have non-trivial 
ground states. We investigate the ground state of quasi-one-dimensional 
(Q1D) Heisenberg antiferromagnets (HAF's) which consist of coupled one- 
dimensional (ID) spin chains with bond alternation on a square lattice. An 
isolated ID uniform spin chain has no long-range order and a striking phe- 
nomenon associated with the Haldane gap [1] is known, namely, chains with 
integer spins have a finite excitation gap over their ground states, whereas 
those with half-odd-integer spins do not. On the other hand, HAF's on a spa- 
tially isotropic square lattice have the long-range Neel order in the ground 
state [2]. What kind of ground states do the models have that lie in the 
intermediate region between genuine ID systems and two-dimensional (2D) 
systems? This is our question. 

The ground state of a ID bond-alternated HAF's with spin magnitude S = 
1 has been investigated extensively [3]. There are two gapped ground states, 
the Haldane phase and the dimer phase, between which a quantum phase 
transition of the Gaussian universality class occurs. We study the ground 
state of coupled 5=1 bond-alternated spin chains on a square lattice. There 
are two ways of coupling such chains, whether we place the stronger bonds on 
parallel positions between the neighboring chains or place them in a zig-zag 
way. The arrangements of bonds are shown in Fig. 1.1. Hereafter we will refer 
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to the former lattice as the square lattice with columnar dimerization and 
the latter as that with the staggered dimerization. The Hamiltonian for the 
model with columnar dimerization is written as follows. 

TL = J ^ (S21-1J • S-2ij + aS 2 ij • S 2 i+ij) + Sij • Sj^+i (1.1) 

The spin operator Sjj has magnitude |S| = 1 and i,j denote the points 
on a square lattice. We consider only the nearest neighbor coupling. We set 
the stronger intrachain coupling unity and the antiferromagnetic interchain 
coupling J' positive. The strength of bond alternation is parametrized by a 
(0 < a < 1). The Hamiltonian of the model with staggered dimerization is 
similarly defined. We set the x axis parallel to the chains. 
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Fig. 1.1. Coupled dimerized chains on a square lattice with (a) columnar dimer- 
ization and (b) staggered dimerization 

There are some previous works on coupled uniform chains [4-6] and 
weakly coupled bond-alternated chains on a square lattice with columnar 
dimerization [6]. It is known that the models have gapped ground states if 
the interchain coupling is weak enough and quantum phase transitions occur 
by tuning the interchain coupling. In our study, by using the continuous- 
time quantum Monte Carlo loop algorithm with the subspin-symmetrization 
technique [7] , we determine the phase boundaries accurately and discuss the 
universality class of the quantum phase transition. As our method is nu- 
merically exact and non-perturbative, we can derive the ground-state phase 
diagram parametrized by a and J' over the whole region. 

1.2 Simulations and Finite-Size Scaling 
1.2.1 Details of Simulations 

Let us denote the size of the lattice simulated by L x , L y and the temperature 
T = 1/ [3. By the Suzuki- Trotter decomposition, we map the original quantum 
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system on L x x L y lattice to a classical system on L x x L y x f3 spacetime by 
adding the imaginary-time axis (denoted by r) with the length (3. The aspect 
ratios L x /L yi L x //3, are fixed. Fixing the ratios between the spatial size and 
the imaginary-time size is based on the assumption of Lorentz invariance [8] , 
namely, the dynamical critical exponent be equal to unity. The sizes of the 
simulated systems and inverse temperatures are up to L x x L y ~ 10 3 sites 
and (3 ~ 10 2 , with 10 3 Monte Carlo steps used for thermalization and 10 4 for 
measurement. The latter are cut into 20 bins from which we obtain averages 
and variances as estimates of observables and their statistical errors. 

We calculate the staggered susceptibility xMi correlation lengths £ x , £ y 
and the excitation gap A. Correlation lengths are calculated from the second 
moment of correlation functions and the gap is obtained as a reciprocal num- 
ber of £ T . By the finite-size scaling (FSS) of these observables, we determine 
critical points and exponents in the thermodynamic limit L Xl L y — > co and 
the ground state limit (3 — » oo. We obtain critical exponents of the staggered 
susceptibility 7 and the correlation length u, by which the critical behavior 
of x(7r) and £<j are described as xM ~ * 7 an d £d ~ where t is the 
distance from the critical point. These are sufficient to give other exponents 
with the scaling relations. 

1.2.2 Determination of critical points and exponents 

We fit the behavior of observables near critical points into the FSS formulae 
which are written as = Lfd{tL x / v ) and xM = L 1 ^ g{tL x / v ), where L is 
the linear system size and fa, g are polynomials. We take terms of the poly- 
nomial up to the second order. Here we describe how the critical point and 
exponents are determined in the ground state of staggeredly coupled bond 
alternating chains with a = 0.1. The raw data of the staggered susceptibility 
Xs are plotted in Fig. 1.2 (a) and its FSS is shown in Fig. 1.2 (b). Data with 
L x = L y = (3 = 16, 24, and 32 are used. As seen in Fig. 1.2, the data near 
the critical point are scaled quite well by choosing J' c , v and 7 as 0.1943(4), 
0.69(2) and 1.4(1), respectively. These exponents coincide with those of 3D 
classical Heisenberg models v = 0.7048(30) and 7 = 1.3873(85) [9] within 
numerical accuracy. It is thus confirmed that the universality class of the 
quantum phase transition of the 2D quantum Heisenberg model belongs to 
that of 3D classical Heisenberg models. The FSS on correlation lengths gives 
consistent results. 

Other critical points are determined in the same way to yield the phase 
boundary over the whole parameter region. 

1.3 Results and Discussions 

In Fig. 1.3, we present the ground-state phase diagrams for columnar and 
staggered dimerization. First of all, both of them arc quite similar with each 
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Fig. 1.2. (a) Raw data plot of the staggered susceptibility and (b) its finite-size 
scaling plot, with L denoting the system size simulated 
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Haldane (y) 


AF 


dimer \ 










\ Hakiiinc (x)" 



Haldane (y) 






AF 








\^ Haldane (x)' 



0.4 0.6 O.S 



Fig. 1.3. Ground-state phase diagrams of the model on a square lattice with (a) 
columnar dimerization and (b) staggered dimerization 



other for large a. Actually these two models are identical on the line with 
a = 1. On this line, there are three points which have been investigated in 
detail so far. The 2D isotropic HAF at (a, J') = (1, 1) has a gapless gound 
state with finite staggered magnetization. The AF phase, which includes the 
isotropic point, occupies a large area in the phase diagram. On the other 
hand, the system consists of decoupled Haldane chains parallel to the x (y) 
axis in the J' = (J' — > oo) limit. The finite excitation gap (Haldane gap) 
observed at J' = and oo survives even at finite J'. We refer to these two 
gapped phases as Haldane (x) and Haldane (y) phases, respectively. 

Now we take a more detailed look on the models with columnar dimeriza- 
tion. At (a, J') = (0, 0), the system consists of decoupled dimers and therefore 
has a spin-gapped gound state (dimer phase). This phase also extends to fi- 
nite a and J'. Most striking feature of the phase diagram shown in Fig. 1.3(a) 
is that the dimer phase around (a, J 1 ) = (0, 0) and the Haldane (y) phase at 
large J' are actually the identical phase. Since the boundary of the AF phase 
does not touch the a = line, there is no critical point between the dimer and 
Haldane (y) phases. Especially it should be pointed out that the a = line 
corresponds to the two-leg ladders, which always has a spin-gapped gound 
state irrespectively of the strength of rung coupling [10]. 

Furthermore, the Haldane (x) and Haldane (y) phases are also shown to 
be identical by considering the bond alternation in the y direction [11]. Thus 
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in the ground states of Q1D models, all of the gapped phases are identical. 
On the other hand, in a strictly ID chain, the dimer phase and the Haldanc 
phase are definitively distinguished in terms of the topological hidden order 
measured by the string-order parameter [12], which is non-zero only in the 
Haldane phase. It should be emphasized that once we have a finite interchain 
coupling the string-order parameter vanishes even in the Haldane (x) phase. 
In this sense the line J' = 0, which represents strictly ID chains, is singular 
in the phase diagram. 

The ground-state phase diagram of the model with staggered dimcriza- 
tion is topologically different from the columnar one. Particularly the AF 
phase extends onto the a = line, as the lattice remains connected two- 
dimensionally even at a — (as long as J' is finite). Thus in the phase 
diagram three spin-gapped phases (Haldane (x), Haldanc (y) and dimer) are 
separated by the AF phase. It is of great interest to pursue phase diagrams in 
the presence of the bond alternation in the y direction both in columnar and 
staggered ways and see the topology of these gapped phases in the extended 
parameter space. 
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